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Abstract 

This paper gives a remark on the Enclosure Method by considering inverse ob- 
stacle scattering problems with a single incident wave whose governing equation is 
given by the Helmholtz equation in two dimensions. It is concerned with the indica- 
tor function in the Enclosure Method. The previous indicator function is essentially 
real-valued since only its absolute value is used. In this paper, another method for 
the use of the indicator function is introduced. The method employs the logarith- 
mic differential with respect to the independent variable of the indicator function 
and yields directly the coordinates of the vertices of the convex hull of unknown 
polygonal sound-hard obstacles or thin ones. The convergence rate of the formulae 
is better than that of the previous indicator function. Some other applications of 
this method are also given. 
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1 Introduction and statements of the main results 



The aim of this paper is to add further new knowledge on the Enclosure Method. The 
Enclosure Method was originally introduced in [5, 4] for inverse boundary value problems 
for elliptic equations. The method aims at extracting information about the location and 
shape of unknown discontinuity embedded in a known reference medium that gives an 
effect on the propagation of the signal, such as an obstacle, inclusion, crack, etc. The 
method can be divided into two versions. One is a version that employs infinitely many 
pairs of input and output data, that is, the Dirichlet-to-Neumann map (or Neumann-to- 
Dirichlet map). Another is a version that employs a single set of input and output data. 
We call this second version the single measurement version of the Enclosure Method. We 
refer the reader to [11, 12, 13] for recent applications of the single measurement version 
of the Enclosure Method. 
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This paper is concerned with the indicator function in the single measurement version 
of the Enclosure Method. In [8], having the single measurement version of the Enclo- 
sure Method, the author considered the reconstruction issue of inverse obstacle scattering 
problems of acoustic wave in two dimensions. The problem is to reconstruct a two di- 
mensional obstacle from the Cauchy data on a circle surrounding the obstacle of the total 
wave field generated by a single incident plane wave with a fixed wave number k > 0. 
The author established an extraction formula of the value of the support function at a 
generic direction which yields information about the convex hull of polygonal sound-hard 
obstacles. However, the indicator function used in [8] is essentially real-valued since only 
its absolute value is used. In this paper, another method for the use of the indicator 
function is introduced. It is shown that the logarithmic differential with respect to the in- 
dependent variable of the indicator function yields directly the coordinates of the vertices 
of the convex hull of unknown polygonal sound-hard obstacles or thin ones. 

Let us describe our main results. First consider a polygonal obstacle denoted by D, 
that is: D C R 2 takes the form D\ U • • • U D m with 1 < m < oo where each Dj is open 
and a polygon; Dj fl D,y = if j ^ f. 

The total wave field u outside obstacle D takes the form u(x; d, k) = e tkx ' d + w(x) with 
k > 0, d E S 1 and satisfies 



where r = \x\ and v denotes the unit outward normal relative to D. The last condition 
above is called the Sommerfeld radiation condition. Some further information about u 
are in order, u belongs to C°°(R 2 \ D) and satisfies u\b £ H l ((R 2 \ D) fl B) for a large 
open disc B containing D. This restricts the singularity of u at the corner of D. The 
boundary condition for du/dv on dD means that D is a sound-hard obstacle and should 
be considered as a weak sense. 

Note that, all the results in this paper k is just positive and there is no other assumption 
on k. 

Let Br be an open disc with radius R centered at a fixed point satisfying D C Br. 
We assume that Br is known. Our data are u = u(-;d,k) and du/dv on OBr for a 
fixed d and k, where v is the unit outward normal relative to Br. Let u and oj 1 - be two 
unit vectors perpendicular to each other. We always choose the orientation of u 1 - and oj 
coincides with e.\ and e 2 and thus oj 1 - is unique. 

Recall the support function of D: ho{oj) = sup^^, x ■ oj. We say that oj is regular with 
respect to D if the set dD fl {x E R 2 \ x ■ oj = ho(^)} consists of only one point. 

Set c t (oj) = tcj + iy/r 2 + k 2 oj 1 - with r > 0. Let v T (x) = e x ' CT ^\ This v satisfies the 
Helmholtz equation in the whole plane. 



Au + k 2 u = inR 2 \D, 



du 
dv 



= ondD, 



r 




Define 




(1.1) 
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This complex valued-function of r is called the indicator function in the single measure- 
ment version of the Enclosure Method. In [8] Ikehata has established the formula: 

lim - log \I(t; uj, d, k)\ = h D (tu), (1.2) 



T 



provided u is regular with respect to D. In this formula one makes use of only the 
absolute value of indicator function I(t;uj, d, k). Thus one needs two regular directions u 
for determining a single vertex of the convex hull of D since formula (1.2) gives only a 
single line on which the vertex lies. Here we present a method for the use of the complex 
values of the indicator function which directly yields the coordinates of a vertex of the 
convex hull of D with indicator functions for a single regular direction u 
Since (uj- 1 -) 1 - = — uj, we have \/t 2 + k 2 uj + iru 1 - = icriu 1 -). This gives 

drVr = i <?\ >? x ' C ^( W± K 

\/t z + k 2 

and thus we have 

d^( drV ^ = ^ T 2 + k 2 { c -( w± ) ■ v + (x ■ c t (uj l ))(c t (uj) ■ u))v T (1.3) 

I\t- u>, d, k) = ■ d T v T - ^(drV T )u^j dS. (1.4) 

Our first result is the following theorem. 

Theorem 1.1. Let ui be regular with respect to D. Let x e dD be the point with 
x ■ uj = ho{oj). There exists a r > such that, for all r > r |/(r; u,d,k)\ > and the 
formula 

I'(t; uj, d, k) 



and 



lim t/ ' L JT\ = h-D(u) +ix -u , (1.5) 



is valid. 



Theorem 1.1 means that, as r — > oo the logarithmic differential of I(r;u,d,k) with 
respect to r converges to x ■ u> + ix$ ■ w -1 . The convergence rate of (1.5) is better than 
that of (1.2). For this see Remark 2.1 in Section 2. 

Remark 1.1. By (1.4) I'(r;u,d,k) can be considered as another indicator function, 
however, it is easy to see that from (1.2) and (1.5) we have 

lim —\og\I\r;u,d,k)\ = h D (uj). 

Thus from this formula we can not obtain any new information about D. 

Here we explain why (1.5) gives further information about D. Let u be regular with 
respect to D. We denote by x(u>) = (x(u)i, x(u) 2 ) the single point in dDn{x e R 2 | x-ou = 
h D (oo)}. Since it holds that 

x ■ (u + iuj^) = {x\ — ix2)(u)i + ioj?), x e R 2 , (1.6) 
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from (1.5) we have 



x(u)i = Re < (u>i + IUJ2) lim 



I'(r; oj, d, k) 



00 \ I(t; oj, d, k) 



and 



x(oj) 2 — Im < (o>i + ioj 2 ) lim 



I'(r; oj, d, k) 



t—,oo y W; d, k) 



;i-7) 



The set /(-D) of all directions which are not regular with respect to D is finite. Let 
1(D) = {(cos 9j, sin 9j) | < Q x < ■ ■ ■ < 6 N < 2tt} and k = 1,2. Define 



x(9) t 



{ x((cos 6, sin 0)) fc , if 6 E [0, 2tt[\{^ \ j — 1, - • ■ , N}, 

i£0 = Oj,j = l,---,N. 



00, 



and extend it as the 27r-periodic function of G R. Since D is polygonal, both x(0)\ and 
x(^) 2 are piece-wise constant and one of which has discontinuity at 6 — 9j in the following 
sense: for each j it holds that x(9j + 0)i 7^ x(6>j — 0)i or + 0) 2 7^ x(0j — 0) 2 . Therefore 
one can expect that computing both x(oj)i and x(oj) 2 for sufficiently many oj via formulae 
(1.7) and (1.8), one can estimate 1(D). This is a new information extracted from the 
indicator function in the Enclosure Method. 

Another implication of Theorem 1.1 is the following idea. Given y e R 2 define 

I(r; y, oj, d, k) = e~ v ' (™+*^T^- ± )/( r; w> rf, fc). 

This corresponds to substitute x i-> e ( :E -J/M™+ l ^' r2 + fc2w± ) instead of v T into (1.1). Since 

l'(r; y, oj, d, k) = -y ■ L + i J—^ uA I(r; y, oj, d, k) + e -H™+*v^W)/' (r . w> d> fc)> 



Vr 2 + A; 2 



it follows from (1.5) that 



v J/ ( r ; y,w,d,k) . ± 

lim - 77 r-rr- = (x — y) ■ (oj + ioj ). 



r^oo I( T ;y,oj,d,k) 
This together with (1.6) yields that 

I'(r;y,u,d, k) 



lim 



I(r;y,u,d,k) 



= Fo -y 



Since rr is the unique point which minimizes the function y 4 |x — y|, one possible 
alternative idea to find x is to consider the minimization problem of the following function 
for a suitable r: 

I'(r;y,uj,d, k) 



y 



I(r;y,u,d,k) 

Since this paper concentrates on only the theoretical issue of the Enclosure Method, we 
leave the numerical implementation of this idea for future research. 
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The result can be extended to a thin obstacle case. Let S be the union of finitely 
many disjoint closed piecewise linear segments Si, S 2 , • • • , S m . Assume that there exists 
a simply connected open set D such that D is a polygon and each Sj consists of sides of 
D. We assume that D C with a i? > 0. We denote by v the unit outward normal on 
dD relative to B R \ D and set u + = v and v~ = — v on S. Given k > and d £ S 1 let 
■u = -u(x), x G R 2 \ S be the solution of the scattering problem 

(A + k 2 )u = in R 2 \ S, 



= on S, 



lim y/r l^r- — ikw] = 0, 

r — >oo \ Qf J 

where w = u — e lkx ' d , u + = m| R 2\75 and u~ = u\d- Note that this is a brief description of 
the problem and for exact one see [8]. Define 



f ( Ou Ov 
Ix(r;uj,d,k) = I —v T - -^u ) dS. 



The following is our second result. 

Theorem 1.2. Let ou be regular with respect to S. // every end points of Si, S 2 , • • • , S m 

satisfies x-u < hs(uj), then there exists a r > such that, for allr > r |/e(t; u,d,k)\ > 
and the formula 

I^(r;uj,d,k) . ± 

lim Tl —r- = h D (u)+ix -uj , (1.9) 

is valid. If there is an end point xq of some Sj such that xq ■ to = hs(u), then, for d that 
is not perpendicular to v on Sj near the point, the same conclusions as above are valid. 

Note that v on Sj n B v (x ) for sufficiently small r\ > becomes a constant vector if 
x is an end point of Sj. 

In [8] under the same assumption as Theorem 1.2, it is shown that 

lim -log|/ E (r;o;,d,A;)| = /i E H- (1-10) 



r 



Thus (1.9) also adds a further knowledge on the use of the indicator function in thin 
obstacle case. 

A brief outline of this paper is as follows. Theorems 1.1 and 1.2 are proved in Sections 2 
and 3, respectively. Both proofs employ some previous computation results done in [8] for 
the proof of (1.2) and (1.10), however, some nontrivial modifications of the computation 
are also required. The idea of using the logarithmic differential of the original indicator 
function developed in this paper can be applied to several other previous applications of 
the Enclosure Method published in [4, 6, 7, 9, 10, 12, 13]. In Section 4 two applications 
of the argument for the proof of Theorem 1.1 are given. In Appendix first for reader's 
convenience we give the proof of Proposition 2.1 which ensures an expansion of the solution 
of the Helmholtz equation at a corner. The proof is focused on some technical part that 
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is different from the case when k — 0. Second a proof of Lemma 2.3 which is important 
for the computation of the expansion of J'(t; oj, d, k) as r — > oo is given. Third a proof 
of some estimates that are needed for the proof of one of two applications in Section 4 is 
given 

Note also that in Sections 2 and 3 we simply write v T — v. 



2 Proof of Theorem 1.1 

Let x denote the single point of the set {x \ x ■ uj = h D {oj)} D 3D. x has to be a vertex 
of Dj for some j. The internal angle of Dj at xq is less than 7r and thus 2n minus the 
internal angle which we denote by satisfies 7r < < 2tt . 

In what follows we denote by Br{xq) the open disc with radius R centered at xq. If 
one chooses a sufficiently small rj > 0, then one can write 

B 2 ri(x ) n (B Rl \D) = {x + r(cos 9 a + sin 9 a 1 ) | < r < 2r], < 9 < 0}, 

B v (x ) n dD = r p u r q u {x } 

where a = cos p u ± + sin p a;, a L = — sin p uj' l + cos p w; — 7r < p < 0; Y p = {x + ra\0 < 
r < T q = {x + r(cos a + sin a L ) | < r < 77}. Note that the orientation of a, a^- 
coincides with that of e 1 ,e 2 . See also Figure 1 of [4]. 

The quantity — p means the angle between two vectors u 1 - and a. p satisfies > 
7r + (— p). Set q = — 2n + p. Then we have — n < q < p < and the expression 

Tp = {x + r(cos pco ± + sin p w) | < r < 77}, 

Tg = {^0 + ?"(cos quj' L + sin gw) | < r < 77}. 

This is the meaning of p and g. 
We set 

u(r, 9) = u(x), x = xq + r(cos 9 a + sin 6* a -1 ). 

The followng proposition describes the behaviour of u(r, 9) as r — )> 0. 

Proposition 2.1 (Proposition 4.2 in [8]). Let rj satisfy rj « l/2k. Then, there exists 
a sequence a±, a 2 , • • • , a m , ■ ■ ■ such that: 
(1) for each s e]0, 2[ 



«(r, 6>) = ^ a m J^ m (kr) cos /i m 6», inH l (B sv (xQ) n \ D)) 

m — l)7r 



m=l 





and J Mm denotes the Bessel function of order ji m given by the formula 



Mm 00 f — I)™ / Z\ 2n 



2/ ^ n!r(n + l+/i m ) V2 
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(2) as m — > oo 

W.o(^(i)1; 

\ y/^m \ r l K J I 

(3) for each I — 1, • • • there exists a positive number C\ such that, for all r e ]0, 77 [ 

|«(r,0)- X> m J Mm (fcr)| < Qr^ 

m=l 

(2.1) 

|u(r,6) - £ « m (-l) m "V Mm (A;r)| < Qr^ 1 , < r < 77. 

m=l 

In [8] the proof is omitted since it can be done along the same line as the proof in the 
case when k = given in [3] . However, there is a technical difference from the case when 
k = and so to make sure and for reader's convenience, in Appendix we give the proof 
which focused on the diffrence. 
Let s = Vt 2 + k 2 + t. We have: 



= H-t)'^=K' + t)- 



Recall that the proof of formula (1.2) is based on the following two facts in [8]. 

• As s — > 00 the complete asymptotic expansion 

e -iV^W e -rfc D («) j(t . w> d> fc) „ £ 62 fc a " |e J I iJ 6 1, (2.2) 

n=2 SM " 

is valid. 

• 3n > 2 a n {e iw " + (-l)V ?/i ™} ^ 0. Thus the quantity 

m* = mm{m > 2 | a m {e ipfMm + {-l) m e iq ^} ^ 0} 

is well-defined, m* depends on k, d, D and u>. 

For the proof of Theorem 1.1 we compute the asymptotic expansion of I'(t; u,d,k) as 
s — > 00. The result is: 

Vr 2 + k^ e - iV ^+^ X0 - U}± e- ThD ^I'(r; u, d, k) 
= -i £ « m {e^- + (-l)- e ^}(- Atm + ^ • Cr (c^)) + O ( — — . 

m=2 15 \o / 

The proof of (2.3) is given in Subsection 2.1. Here we show how to prove Theorem 1.1 
by assuming (2.3). 

Since a m {e ip ^ m + (-l) m e ^™} = for all m with m < m* and (3 = a m *{e ip ^ m * + 
(-1)™* e^™*} ^ 0, from (2.2) we have 

/ K\ Mm* 

^lnn^ t^* e -*/^^W e -rftz>M j( r; ^ ^ fc ) = _ i/3e *§M m * _ ( 2 .4) 
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On the other hand, since 



ix ■ c T (u r ) . ± 

lim . = Xn ■ CO + tXn • CO , 

r^oo Vr 2 + k 2 

it follows from (2.3) that 

/ 7 \ Mm* 

lim r^*e- lv/7Wa;o '" X e- r ' lo(a;) /'(r;a;,d,A;) = -^(xo-w+ixo-w^e^^' . (2.5) 

r >oo V 2 / 

A combination of (2.4) and (2.5) ensures the validity of (1.5). This completes the proof 
of Theorem 1.1. 

Remark 2.1. Since Crioo^) depends on r and thus s, (2.3) is not the complete asymptotic 
expansion. However, using (2.2), (2.3) and the expression 

s k 2 
ic T (co L ) — -(to + ico L ) + —(to - ioo L ), 

one can easily obtain the following expansion: 

= i V a m {e ip ^ + (-1)™^}^ — + O ( — — 

where I = J(t; to, d, k) and I' = J'(t; w, d, fc). Note that we have used n < < 27r. 

The above formula yields the second term of the expansion of the logarithmic differ- 
ential of the indicator function as r — > oo: 

I'(r;uj,d,k) . ± /x m . 

I(r;oo,d,k) = M")+«*o-* - — + y T 2, 

Thus the convergence rate of (1.5) is better than that of (1.2) since (2.4) yields 
- log | J(r; w, d, fc) | = M") - f^^L + O (-) . 

T T \T / 

From this view point one can say that (1.5) is an improvement of (1.2). Note also that 
both formulae show that the accuracy of the approximation depends on the size of /i m * = 
(m* - l)7r/e. 

2.1 Proof of (2.3) 

Integration by parts gives 

d 



f o 

I'(t; to, d,k) = / u—d T vdS. 

JdD ov 
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Localizing this integral at xq, we have, modulo exponentially decaying as r — > oo 



~ e 



iV^+^x - w - L e -rfc D H /* u Q rV dS + e -iy/i"+Wxo+r L e -Th D (u) f u d T vdS (2.6) 
Jv v ov Jr a ov 



= I p (r) + I q (r). 

We have:// = sin pus 1 - — cos and x = Xq + r(cos pa;- 1 + sin pco>) on T p ; v 
— sin go;- 1 + cos qu and x = Xo + r(cos gw 1 + sin qui) on T g . 
From those we have: 
on r p 

Criuj- 1 ') ■ v = r sin p + ix/ r 2 + k 2 cos p, 



c T (u;) • 1/ = — (r cos p — iV r 2 + A; 2 sin p), 
x • c^w -1 ) = Xq ■ c T (co ± ) + r(r cos p — Wt 2 + k 2 sin p), 



r/i£)(u;) i%/ r 2 +fc 2 2!o-w x „Tr sinp ix/r^+fc^Y cos p. 



on r 



c^w 1 ) • 1/ = — (r sin g + i\J r 2 + A; 2 cos g), 
c T (w) • i/ = r cos g — iy/ r 2 + A; 2 sin g, 
x • c^w -1 ) = x ■ c T (co ± ) + r(r cos g — i\Ar 2 + k 2 sin g), 



p rh]j (uj) i\f r 2 +k 2 xq-ui 1 - rr sin q i\/ t 2 +k 2 r cos q 



(2.7) 



(2.8) 



It follows from (1.3) and (2.7) that 

Vr 2 + k 2 I p (r) 

= i |(r sin p + ^Vt 2 + A; 2 cos p) — x • c t (co ± )(t cos p — r 2 + A; 2 sin p)| 

x f V u(r,0)e Trsinp e iVTTT ^ rcosp dr ( 2 - 9 ) 

JO 

-i(r cos p - iVr 2 + A; 2 sin p) 2 /" ru(r, 0)e rr sinj, e iWr cos p dr. 

Jo 
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It follows also from (1.3) and (2.8) that 

V^T¥l q (r) 

= —i |(r sin q + i\/ r 2 + k 2 cos q) — x ■ c r (w L )(r cos q — iy/r 2 + k 2 sin g)| 

r> ■ -act ( 2 - 10 ) 

x / ii(r,e)e Trsm V^ rms «(ir 
Jo 

+i(rcos g-«Vr 2 + fc 2 sin g) 2 /" ru(r, 0)e rr sin V v/?Wr cos 9 rfr. 

JO 

Here we make use of (2.1). Since ot\ J (k\x — x \) satisfies the Helmholtz equation in 
the whole plane, the indicator function for u coincides with that for u — aiJ (k\x — | ). 
Thus one may assume, in advance a± = in the computation of the integrals in (2.9) and 
(2.10). 

Since p and q satisfies — n < q < p < 0, we have sin p < and sin q < 0. This gives, 
for 6 = p,q and /i > 

[ V r^e Trs ' m6 dr = 0(t-^ +1) ). (2.11) 
Jo 



Set 



I^r, 6) = f V J^kr)e TT sin V v/ " Wr cos e dr, 

J 



K^t,0) = I rJ^kr)e TrsAne e iV7ITWIrcose dr. 

J 

It follows from (2.1) and (2.11) that 



n 



(2.12) 



r u {r,0)e Trsinp e i ^ TT ^ rcosp dr = £ a m /^,(r,p) + 0(r-^+ 1+1 )), 

J ° m=2 

[ V ra(r,0)e rrsin V V7Wrcosp rfr = £ a m AT^ m (r,p) + 0(r-^+ 1+2 )), 

"' m=2 



/\(r,0)e^ sin V^ Wrcos ^r = £ a m (-l) m -% m (r, q) + 0(r~^+ 1+1 )), 

[ V ru(r, B)e Tr ^ V^Wr cos <?rf r = V ^(-l)" 1-1 ^ (r, q) + 0(r-^+ 1+2 )). 
Jo ^2 

Substituting these into (2.9) and (2.10), we obtain 
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Vr 2 + k 2 I p (r) 

= % |(r sin p + iy/ r 2 + k 2 cos p) — • c T (w L )(r cos p — i\/ r 2 + k 2 sin p)| 

n 

H a mI„ m {T,p) 



n 

X 

m=2 



(2.13) 



-i(r cos p - iVr 2 + /c 2 sin p) 2 ^ a m K^ m (T,p) + 0(r ^ n+1 ) 



m=2 

and 

Vt 2 + k 2 I q (r) 



—i |(r sin g + i\J t 2 + k 2 cos g) — rr • c r (w _L )(r cos g — i\/r 2 + k 2 sin g) j 

xEM-l) m %(M) 



(2.14) 



ra=2 



+^(r cos g - i\J r 2 + /c 2 sin g) 2 £ a m (-l) m -% m (r, g) + 0(t-^). 



m=2 



In what follows we choose an arbitrary n and fix. Let /x = fi m with < m < n. We 
have already established that 

Proposition 2.2([8]). As r — > oo we have 

(r cos - iVr 2 + A; 2 sin 0)/ m (t, 0) = + 0(s~°°). (2.15) 

The main problem is to compute the asymptotic expansion of the quantity (r cos 6 — 



Wr 2 + k 2 sin 9) 2 K^r,9). 
We prove 

Proposition 2.3. As s — > oo we have 



(r cos 9 - iVr 2 + k 2 sin 6) 2 K^t, 6) 

= <(i - o 2 {i + c + mi - c)}(i - cr 3 — - — + o(a-°°), 

w/iere ( = (k/s) 2 e 2%e . 

For the proof we prepare the following three technical lemmas. 
The first one can be proven along the same line as Lemma 3.1 in [8]. 



(2.16) 
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Lemma 2.1. For each I' — 0, 1, • • • we have 



Y , / ( -) ^ r 2.+M+l e -sin ^.v^TFrcos ^ 

^ oJ !r(i + j+/i) Jo (2-17) 

1 



+0 



r 2+2(/'+l)+/i 



The second is nothing but Lemma 3.2 in [8]. 
Lemma 2.2. Let —n < 9 < 0. For eac/i / = 0, 1, • • • as s — > oo i/ie formula 

[°° r^p-rrsinO iVr*+)?r cos 9 j _ ST^ L^ n (6) ( 1 \ Icy q\ 

is m/id where 

L a , n (6) = ie w e^ + ^2° + \-k 2 e™T V{a + n+l \ 

The third corresponds to Lemma 3.3 in [8], however, the proof needs a careful modi- 
fication of that of Lemma 3.3. See Appendix for the proof. 

Lemma 2.3. 

(-irr(n + 2 + , i2 + ,„,) = 
».Wr(l + n.+fc) w 1 ' 

Proo/ of Proposition 2.3. A combination of (2.17) and (2.18) for a = 2j + /i + 1 gives 

^j!r(l +j +/i) \2J S 2j+ M +2n+2 ^ ^2j+ M +2(/+l)+2y J 

+0 ( s 2+2(2'+l)+„) (2.20) 



i !r (! + J + /i) S 2(n +J -) +M+ 2 



s 2(Z+l)+Ai+2 y \ vs 2(2' + l)+ A t+2 

Now let I = I'. Then (2.20) becomes 
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(-1)™ 2 fk\ 2n2+fl 
= £ LL^r(l + n 2 + ,) UJ ) X (2.21) 



+0 



1 



s 2(z+i)+^t+2 y • 



Write 



(-IT ^V W2+ ^ ^ 



ni+n 2 



(_l)n 2 

"n 1+ „^!r(l + n 2 + /i) 1,2 J 

x ^ c ^ c »(g+f )(2» 2 +M+l)o2n 2 + M +2/ k 2 e 2i6\ni r ( 2n 2 + ^ + ^1 + 2) 

ni! 

= 2 2 ie ie i(e+f)(M+i) fc M/_ fc 2 2iONn (-l)" 2 r(n + 2 + n 2 + fi) 

ni &=n n 1 !n 2 !r(l + n 2 + /i) 



(2.22) 



From (2.19), (2.21) and (2.22) one gets 

Z (V,2 a 2id 

E 

n=0 

From this one gets 

r cos 9 - iVr 2 + k 2 sin #) 2 i^(r, 0) 
_ (se~ w ) 2 



2 2 



■(i-C)' 



x /2 2 < e * e ^5)^^ T (^)> + l)(n+l+^) / 1 



= e—(l-C) 2 

x L-e^^ £ (^)"(n + l)(n + l + ^) / I 

j ^ s 2n+^ ^ s 2(/+l)+/i 

r i (k 2 e 2ie ) n (n + l)(n + l + pi) i ^ ( 1 \\ ieW+^W 

~2 



(2.23) 



- 'I 1 U ^ s 2n ^ U \ s 2(l+l) J ( sM 
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Let |C| < 1. Since 



E O + 1) = (1 " Cr 2 , E O + I) 2 = (1 + 0(1 " 0" 



, V"- ~r *-) -^Tsyv 1- sy 3 ) 

n=0 n=0 

we have 

oo 

£ C> + l)(n + 1 + ii) = {1 + C + Ml " C)}(1 " C)- 3 

n=0 

and for each fixed I 

E C(n + l)(n + 1 + /i) = {1 + C + Ml - OKI - 0" s + 0(|C| ,+1 ). (2.24) 

n=0 

Substituting (2.24) with C = (k/s) 2 e 2ie into (2.23), we have 



(t cos 9 - iVr 2 + k 2 sin 9) 2 K^r, 9) 

= i(i - 2 {i + C + Mi - OKI - cr 3 ^A + o (j^) . 

Since / can be arbitrary large whenever n is fixed, /x = /i m and < m < n, we obtain 

(2.16). 

□ 

We continue the proof of (2.3). One can write 

r sinfl + ^v^TP cos^ = ^^|l+ Q e 2i * 

and 

T C os - ^v^TP sin 9 = ^ j 1 ~ (j) e ™ 
From these and (2.15) one gets 

|(r sin + iVr 2 + A; 2 cos 0) — rr • c^w -1 )(r cos - iVr 2 + A; 2 sin 0)} 7 M (r, 0) 

(2.25) 

, 7 >i(64-fW,u V y 

= {i(l + 0(1 " 0" 1 - ^0 • CM)} + 0(s~°°). 

Since 

i{\ + o(i - 0" 1 - <(i - 2 {i + C + Mi - 0}(i - 0~ 3 = -in, 

it follows from (2.13), (2.14), (2.16) and (2.25) that 

Vt 2 + k 2 I p {r) = V a m (ifi m + x ■ c^ur 1 )) h O — — 

and 



v^TF/,(r) = E a m (-l) m (Vi m + • c^)) + O 

m=2 S V 



Now from these and (2.6) we obtain (2.3). 
□ 
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3 Proof of Theorem 1.2 

First consider the case when every end points of Si, S 2 , • • •, S m satisfies x ■ u < hs(u). 
Then x G S with xq-oj = h^(uj) should be a vertex of D and a point where two segments 
in some S^ meet. We take the same polar coordinates as those of Section 2. 
Integration by parts gives 

d 



[u]—d T vdS, 



where [u] = u + \qd — u \qd- Localizing this integral at xq, we have, modulo exponentially 
decaying as r — > oo 

^-iVWxo^g-r^H / yla T D(iS + e- iV?TF ^ ui e-*( u »/ [«]^-9 T udS' (3.1) 

•/r p ai/ Vr, az/ 

= J p (r) + J 9 (r). 

It follows from (1.3) and (2.7) that 

Vt 2 + k 2 I p (r) = i j(r sin p + r 2 + A; 2 cos p) — rr ■ c T (w" L )(r cos p — r 2 + A; 2 sin p)} 

x f\u + (r,0) -u-{r,2n))e Trsinp e iV7TT ^ rcosp dr 
Jo 

-i(r cos p - iVr 2 + A; 2 sin p) 2 f r{u+{r, 0) - u~ (r, 2n))e rr sin? e ,Wr cos p dr. 

JO 

(3.2) 

It follows also from (1.3) and (2.8) that 
Vt 2 + k 2 I q (r) = — i |(r sin g + r 2 + A; 2 cos g) — Xo • c T (w _L )(r cos g — iV r 2 + A; 2 sin g) j 

x /% + (r, 0) - u-(r, 0))e rr sin V v/ " Wr cos q dr 
Jo 

+i{r cos g - iVr 2 + A; 2 sin g) 2 /" r (w + (r, 0) - u" (r, 0))e rr sm^VTWr- cos 

JO 

(3.3) 

By Proposition 4.4 in [8], we have 

oo 

u +(r, 0) = ^ a+ J M + (Ax) cos < r < 77, < 9 < 0, 

m=l 

00 

«"(r, 0) = ]T a" J Mm (Ax) cos ^(9 - 0), < r < 77, < 9 < 2n, 

m=l 

where //+ — (m — l)7r/0 and /x~ — (m — l)ir/(2ir — 0). For the precise meaning of this 
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expansion see [8]. From these and /x+ < ji m we have 



i 



[u] p = u+(r,0) - u-(r,2n) = £{J,+ (A;r)a+ + J , m (kr)a m {-l) m } + O(r^), 



m=l 
I 



[u] q = u + (r, 0) - u-(r, &) = -J2 {J^{kr)at{-\r + J, m (kr)a m } + O(r^). 

m=l 

It follows from (3.4) and (2.11) that 

f\u + (r, 0) - «"(r, 27r))e Tr sin V v/7Wr cos p dr 

JO 

n 

= E KW^p) + (r,p)} + 0(r-^ + +1 +1 )), 

m=l 

/" r(« + (r, 0) - u"(r, 27r))e rr sin P e i ^ T +^ r cos *>dr 

JO 

= E fop) + ^(-1)^,- (^P)I + 0(r~^ +1 +2 )), 

m=l 

/% + (r, 6) - «"(r, 0))e rr sin V v/7Wr cos grfr 

JO 

ra 

= " E{«m(-l) m ^(^?) + «m^(^?)} + 0(r-^ + i +1 )), 
m=l 

fr(u + (r, 6) - «"(r, 6))e rr ^ v ^Wr cos grfr 

JO 



- E {^(-l) m ^(^ ?) + ^-Jt, q)} + 0(r-« +1 



+2h 



m=l 



Substituting these into (3.2) and (3.3), we obtain 



Vt* + WI p {t) 

= i j(r sin p + i^r 2 + /c 2 cos p) — x ■ c t (co j ~)(t cos p — iy/r 2 + k 2 sin p)} 

n 

E{«^(^p) + <(- 1 )%^(^p)} 



n 

X 

m=l 



z(t cos p - i\/r 2 + /c 2 sin p) 2 £ (r,p) + a" (-l) m ^ m (r,p)} + 0(r 



m=l 
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and 

Vr 2 + fc%(r) 

= sin g + r 2 + /c 2 cos g) — x ■ c t {u l )(t cos q — iV r 2 + k 2 sin g 



x E {a+(-l)%+ (r, g) + a"/^(r, g)} 

m=l 

n 

+, (r cos g _ ?V ^TF sin g) 2 E {^(-1)"^+ (r, g) + a~K^- (r, g)} + 0(r 

m=l 

It follows from (3.5), (3.6), (2.16) and (2.25) that 

Vr 2 + fc 2 / p (r) = E + ^0 • ^(W- 1 -)) "+ 

m=l sMm 

n p i(P+^)lJ-mU^m ( 1 \ 

+ E « m (-l) m (^™ + • c r (^)) + O "H 

and 

v^ + F/ ff (r) = E a+(-l) m Mi + • c^)) -+ 

m=l ° 

n e *(g+f )Cmpm / 1 \ 

+ E a m(^m + ^0 • C r (^ ± )) + O — — . 

m=l « Mm Vs^+i/ 

Now it follows from this and (3.1) that 

= ~i E «-{ eWm + (-l) m e^-}(-^ + is • cV)) 

m=l 



-< E «-{(-l) m e^ + e^}(-/i m + *:r • c r (c^)) — + O 

m=l S ^ m VS-™+ 

Since we have the following cancellation ([8]) 
(3.7) becomes 

= -» E «m{ eiPMm + (- 1 ) m e^ m }(-A i m + • c T (w ± )) + + O " 



m=l 
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By the way, we have already known in [8] that 

= -iV a+{e^ + (-l) m e <w *} - + O 

m=l ^ m 

and there exists a m > 2 such that a+{e w ™ + (— 1)"V 9M ™} 7^ 0. Having these together 
with (3.9), hereafter we take the same course as the obstacle case. 

The case when there is an end point of some Ej such that x ■ u — h^{oj) corresponds 
to the case when p — q. We omit its description. 
□ 

Remark 3.1. From (3.9) and (3.10) we see that the field u~(r, 9) never affects the asymp- 
totic behaviour of I^(t;oj, d, k) and /s (r; u>, d, k) as r — > 00 modulo rapidly decreasing. 
The key point is cancellation (3.8). 



4 Some other applications 

In this section we present some implications of the argument done for the proof of Theo- 
rems 1.1. 



4.1 From the far- field pattern of the scattered wave for a single 
incident plane wave. 

Let u = e lkx ' d + w be the same as that of Theorem 1.1. It is well known that w has the 
asymptotic expansion as r — > 00 uniformly with respect to ip e S 1 : 

W {r V ) = e —F{^k) + 0[-^). 

The coefficient F((p; d, k) is called the far-field pattern of the scattered wave w at direction 

ip. 

In this subsection we present a direct formula that extracts the coordinates of the 
vertices of the convex hull of unknown polygonal sound hard obstacles D = D x U ■ • • D m 
from the far field pattern for fixed d and k. 

We identify tp = (pi,<p 2 ) with the complex number given by tpi + iip2 and denote it 
by the same symbol (p. Given N = 1, • • •, r > 0, 00 G S* 1 and k > define the function 
<7at( ' ; T , k, u>) on S* 1 by the formula 



1 ^ f ikip 



Then we have 



d T g N (tp ] T,k,u) = -— 7 === £ 



ikip 



2nV^+¥ (r + V^+¥) 



CO 
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In this subsection Br denotes the open disc centered at the origin of the coordinates with 
radius R and we assume that D C Br. 

Theorem 4.1. Let u> be regular with respect to D. Let f3 be the unique positive solution 
of the equation 

2 

-s + log s = 0. 

e 

Let (3 satisfy < j3 < (3 . Let {t(N)} n=1> ... be an arbitrary sequence of positive numbers 
satisfying, as N — > oo 

T(J v ) = ^ + 0(1 ). 

Then the formula 



F(-ip; d, k)d T g N (ip; t(N), k, u)dS(ip) 

lim = hr>(uj) + ixo ■ cu" 1 , (4.1) 

/ F(-ip; d, k)g N (ip; t(N), k, w)dS(<p) 
J s 



N — >oo 



is valid. 



In [10] we have shown that, under the same choice of t(N) and u> being regular with 
respec to D 



lim — — - log 

TV^oo T (N) B 



S 1 



F(-<p; d, k)g N (ip; t(N), k, uj)dS(ip) 



h D (u). 



Thus Theorem 4.1 corresponds to Theorem 1.1. 

Let us describe the proof of Theorem 4.1. The starting point is the following identity 
which is a consequence of the formula (2.9) in [2]: 

V87rk 



in/4 I F(-(p;d,k)g N ((p;r,k,u)dS(cp) 
e 1 Js 1 

(4-2) 

= J(t; u, d, k )+ J 9B | ~^( v 9n - v r) ~ q^( v 9n - v t)u I dS 

where 

v g N (y) = e ty " p g N (<f,T,k,uj)dS(<f). 



s 1 



By Theorem 2.1 in [10] we know that the second term in the right-hand side of (4.2) has 
the bound e- T ^ R 0(N-°°) as N — ► oo. Since e- T (*H = 0(e T W R ), it follows from 
(2.4) and (4.2) that 



lim T (iV)^ e-iv^W-o'^g-r f F(-<p;d,k)g N (<p;T(N),k,uj)dS(<p) 



(kY 



2 J 

(4.3) 
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Similarly to (4.2) we have 

pin/4 



F(-(p; d, k)d T g N (ip; r, k, u)dS((p) 



I'(t; u, d, k) + j q^( V9 t9n ~ d rV T ) - -^{v dT9N - d T v T )u \ dS 



v d T9N {y) = / , e w '' p d T g N (ip;T,uj,k,d)dS((p). 



(4.4) 



where 



The following lemma corresponds to Theorem 2.1 in [10] and see Appendix for the 
proof. 

Lemma 4.1. We have, as N — > oo 



e RT W sup 
\y\<R 



I e ik ^d T g N (cp; r(N), k, u)dS(<p) - d T v T (y)\ 
J s 



t=t(N) 



+e Rr(N ^ sup 
\y\<R 



V | jf si e ik ^d T g N { V] t(N), k, ou)dS(<p) - d T v T (y)\ T=r{N) 



= 0{N-°°). 

Now from Lemma 4.1, (4.3) and (2.5) we obtain 
lim rf^^e-'V^^^e-W^H^? / F(-<p;d,k)d T g N ((p;T(N),k,u))dS((p) 

/ U \ 



■> f 1 n 



= i/3(x ■ u + ix ■ u )e l 2 

y a i 

From this together with (4.4) yields (4.1). 

4.2 From the Cauchy data of the scattered wave for a single 
point source 

Let y E R 2 \ D. Let E = E D (x,y) be the unique solution of the scattering problem: 

(A + k 2 )E = in R 2 \ D, 

'dE 



where 



lim ( ikE I = 0, 

r — >oo \ Q r I 



%(x,y)= l -Hl ) 1 \k\x-y\ 
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and denotes the Hankel function of the first kind [14]. 

The total wave outside D exerted by the point source located at y is given by the 
formula: 

$ D (x, y) = $ (x, y) + E D (x, y),xeR 2 \D. 

In this subsection we consider the following problem. 
Inverse Problem. Let Ri > R. We denote by B R and B Rl the open discs centered at a 
common point with radius R and Ri, respectively. Assume that D C B R . Fix k > and 
y G dB Rl . Extract information about the location and shape of D from $ D (x,y) given 
at all x G dB R . 

Define 

J(r;u,y,k) = ^-§ D (x,y) -v t (x;uj) - ^-v T (x;u) • ® D (x,y)^ dS(x). 
Then we have 

j'{T;u,y,k) = (j^§ D (x,y) ■ d T v T - ^-(d T v T ) ■ $ D (x,y)J dS(x). 

Note that (d/dv)<& D (x,y) for x e dB R can be computed from & D (x,y) given at all 
x G dB R by solving an exterior Dirichlet problem for the Helmholtz equation. See [12] 
for this point. This remark applies also to [d/dv)u on dB R in Theorems 1.1 and 1.2. 

A combination of the proof of Theorem 1.2 in [12] and the same argument as done in 
the proof of Theorem 1.1 yields the following formula. 

Theorem 4.2. Assume that 

diamD < dist(D, dB Rl ). (4.5) 

Let u be regular with respect to D. Let x G 3D be the point with xq ■ u = Hd{oj). Then, 
there exists a r > such that, for all r > r \J{r; u>, d, k)\ > and the formula 

J'(r;u,y,k) . ± 

hm — — = h D (uj) +ix -uj , 

^-^oo J( T ;u,y,k) 

is valid. 

Note that it is an open problem whether one can drop condition (4.5). For more 
information about this see [12]. 
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5 Appendix 

5.1 Proof of Proposition 2.1 

Proof. For each m define 

u m {r, 0) = / u(r, 6) cos // m 0d0, (r, 6) G ]0, 2 V [ x [0, 0]. 
Jo 
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Then we see that u rn satisfies the equation 

-(ryj + (Ek _ k 2 r)y = in]0, 2 V [ 
r 

and thus this yields that there exist numbers a m , f3 m such that 

u m (r) = a m J^ m (kr) + (3 m Y llm (kr), 

where Y^ m (kr) denotes the Bessel function of the second kind ([14]). Then a similar 
argument done in [3] and the behaviour of J^ m (kr) and Y^ m {kr) as r — > ([14]), one 
concludes that f3 m = for all m > 1. Substituting u m (r) = a m J^ m (kr) into the inequality 



/ r\u' m (r)\ 2 dr < \\Vu\\ 
Jo 



r-277 

' 1 2 



i 2 (B 2 ^(^o)n(s fl \D)) 



we obtain 

/•2n 
2 ' 

a 



/ r|{J^(fcr)}'| 2 dr < C (Al) 

where C is a positive constant independent of m. 

Thus the problem is to estimate the integral in the left-hand side of (A.l) from below. 
For the purpose we make use of the following formula which can be checked directly: 



{J»Jkr)}' = ^J»Jkr)-U» m+1 (kr). 



From this we have 



r\{J»Jkr)}'\ 2 = ^\J,Jkr)\ 2 - 2» m kJ»Jkr)J» m+1 (kr) + k 2 r\ J, m+1 (kr)\ 2 
and thus 

r2r) f2f] \ r2r\ 

\ r\{J, m {kr)Y\ 2 dr>^ m \J^kr)\ 2 dr - 2^ m k J,Jk)J, m+1 (kr)dr. (A2) 

JO JO V Jo 

By formula (37) on p. 338 in [1] and a change of independent variable we have 

r2rj r2r]k 

k J^ m (kr)J^ m+1 (kr)dr = J^ m (r)J^ m+1 (r)dr 



(A3) 

oo 

= EIW(¥)I 2 

n=0 

Furthermore from (A. 3) and the recurrence relation of the Bessel functions 

HmJ» m (kr) = fcr(J Mm _i(A;r) + J^ m+1 (kr)) 
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which can be checked directly we have 

r 2r i 1 k f 2r i k r 2r i 

/ -\J»Jkr)\ 2 dr = — J» m -i(kr)J flm (kr)dr + — / J^(kr) J^ +1 (kr)d7 
Jo r ii m Jo fi m Jo 

j r2i)k | r2r)k 

li m JO li m JO 



i oo i oo 

= — E l^ + n(2^)| 2 + — E \J» m+ n + l(2 V k)\ 2 

1 2 00 

= — | J Mm (2^)| 2 + — \J, m+n+ i(2vk)\ 2 . 

This together with (A. 2) and (A.3) yields 

/ r|{J, m (A;r)}frfr> / i m |J, m (2r ? A;)| 2 . (A4) 

Since 

JMm( " ^"ra + Ml + i r(i + ^ + n) w | 

and 

r(i + /x m ) = i < i_ 

r(l + ^ m + n) n^ =1 (j+/i m ) " n!' 

it holds that 

Thus choosing r\k in such a way that e r ' k — 1 < 1/2, that is, r^/c < log(3/2), we obtain 

1 (rjky 



2r(i + /i m )- 

This together with (A. 4) gives the following estimate: 

f 2 V 1 

l r\{J, m (kr)}fdr > - ^ ^ )2 (r)k) 2 ^. (A5) 

Then (A.l) and (A. 5) give the estimate in (2). The remaining parts of the statements are 

consequences of (2) and the completeness of (cos /x m #)^ =1 in L 2 (]0, 0[). 

□ 

5.2 Proof of Lemma 2.3 

Since 



r(n + 2 + 1 + = m+l(j + 1 + /i)}r(i + i + 
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one can rewrite 



^ __ n m!n 2 !r(l + n 2 + y) 



ni+«2 : 



( — l) n2 



1 / /V \ ™ +1 



n! 



n + 1 f d 



= (-l)"(n + l)(n + l + /i). 

□ 



5.3 Proof of Lemma 4.1 

In [10] we have already known that 

v gN (y)-vAy) 



m>7V I fc J m>iV I fc J 

(A6) 

where y = (rcos #,rsin 0). Since 



<9 T 



(r =F Vt 2 + A; 2 )wl ™ f (r =F Vi" 2 + fc 2 )u|™ 1 (Vr 2 + k 2 =F r)ZJ 



= m 



fc j Wr 2 + k 2 



Vt 2 + k 2 
from (A. 6) we have 



m j (r — Vt 2 + k 2 )oo 



Vt 2 + k 2 (v dT g N (y) - d T v T (y)) 



imB 



Since 

|^n(fcr)| < 



A;r\ m 1 
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it follows that 



V^T¥\v dTgN (y) - d T v T (y) \<C(N + l)E(r; N + 1), (A7) 
where C > is independent of N and r and 



N 



E( T ;N) = ± f g(l±Vp±g) j e %+VW)/ 2i 

By virtue of the choice of t(N), (3 and the Stirling formula (cf. [14]), we have 

e Rr{N) E(t(N)- N + 1) = C^iV" 00 ). 

This is the key point and see [10] for the detail of the derivation. This together with 
(A. 7) yields the half of the desired estimates. The remaining estimate can be also given 
similarly. 

□ 
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